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Abstract. The Riemann zeta function at integer arguments can be written as an infinite sum of 
certain hypcrgcometric functions and more generally the same can be done with polylogarithms, 
for which several zeta functions are a special case. An analytic continuation formula for these 
hypergeometric functions exists and is used to derive some infinite sums which allow the zeta 
function at integer arguments n to be written as a weighted infinite sum of hypergeometric 
functions at n — 1. The form might be considered to be a shift operator for the Riemann zeta 
function which leads to the curious values C F (0) = J (2) - 1 and ( F (1) = Ei(l) - 7 which 
involve a Bessel function of the first kind and an exponential integral respectively and differ from 
the values £ (0) = -A and £ (1) = 00 given by the usual method of continuation. Interpreting 
these "hypergeometrically continued" values of the zeta constants in terms of reciprocal common 
factor probability we have £ F (0) _1 = 78.15% and C F - 75.88% which contrasts with the 

standard known values for sensible cases like ( (2)~ 1 = S 60.79% and C(3) _1 = 83.19%. 
The combinatorial definitions of the Stirling numbers of the second kind, and the 2-restricted 
Stirling numbers of the second kind are recalled because they appear in the differential equation 
satisfied by the hypergeometric representation of the polylogarithm. The notion of fractal strings 
is related to the (chaotic) Gauss map of the unit interval which arises in the study of continued 
fractions, and another chaotic map is also introduced called the "Harmonic sawtooth" whose 
Mellin transform is the (appropritately scaled) Riemann zeta function. These maps are within 
the family of what might be called "deterministic chaos" . Some number theoretic definitions are 
also recalled. 

Contents 

1. The Zeta Function 2 
1.1. Riemann's Function 2 

2. Number Theory, Continued Fractions, and Fractal Strings 6 

2.1. Fractal Strings and Dynamical Zeta Functions 6 

2.2. The Gauss Map h(x) J7 

2.3. The Harmonic Sawtooth w(x) 12 

2.4. The Prime Numbers 14 

3. Analytic Continuation 15 

3.1. Continuation ol n+ iF n Near Unit Argument ljj 

3.2. The Continuation of Li^(i) and£ F (ri) via Contiguous Functions 15 

4. Appendix 2C 

4.1. The Generalized Hypergeometric Function : p F q 2C 

4.2. Other Special Functions |21 



Email: Stephen. crowleySmavs .uta. edu. 



1 



STEPHEN CROWLEY 

4.3. Notation [22 
References 22 



1. The Zeta Function 

1.1. Riemann's £(t) Function. Riemann's zeta function, named after Bernhard Riemann(1826- 
1866), is denned by 



(1) =T^E^o(2« + ir t V{teC 

= (1 - 2- f ) ry(t) V{teC 



H(t) > 1} 

ft(t) > i} 
n[t) > o} 



where lZ(t) and denote real and imaginary parts of t respectively and n (t) is the Dirichlet eta 
function, also known as the alternating zeta function, named after Johann Dirichlet (1805- 1859) 

, 2) V(t) =En=i^ Zl V{te€:K(t)>0} 

= r^y / X s - 1 -^d.T V{i G <D : K(t) > 0} 

where the integral is a Mellin transform of (e x + 1) _ . The function is analytic and uniformly 
convergent when *H(£) > 1 or > when using the eta function form. The only singularity of 
is at t = 1 where it becomes the divergent harmonic series. The reflection functional equation 
pi 13.151] which relates CO) to C(l - t) is given by 

(3) ctoTr-^-'rWcosfy) =C(i-*) 

The interpretation of zeta in terms of frequentist probability is that given n integers chosen at 
random, the probability that no common factor will divide them all is £ (n) . In other words, 
given an array i of n random intgers, C ( n ) is the probabability that gcd (ii, 12, ■ ■ ■ , i n ) = 1 where 
gcdis the greatest common denominator function. So for example, the probability that a pair of 
randomly chosen integers is coprime is C (2) 1 = — 60.79%, and the probability that a triplet 
of randomly chosen integers is relatively prime is C (3) _1 = 83.19%. [37] 051 13.1] 1.4] 

1.1.1. The Generalized Hurwitz Zeta Function a). A more general function which includes Rie- 
mann's Zeta function was defined by A. Hurwitz. 

00 

(4) ({t, a) = £> + a)-< 

n=Q 

Notice that the summation starts at n — whereas Riemann's starts at n = 1. It is apparent that 
Q(t) is a special case of Q(t, a) where 

00 00 

(5) c( t ) = ]Tn- t = c(M) = £("'+ i r f 



mi pi 13.11] 
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1.1.2. Hypergeometric Representations of the Lerch Transcendent: <&(z,t,v). The Lerch transcen- 
dent <&(z,t, v) 10, 1.11] is a further generalization of the Hurtwitz zeta function 



(6) *(*>*,«) = T,n=0 

valid V|z| < 1 or {\z\ = 1 : K(t) > 1} which is related to C(f, v) and ((t) by 

*(l,t,w) =C(*,w) 

(7) $(l,t,l) =C(*) 

*(l,t,l/2) =C(t,l/2) =(2*-l)C(t) 

When i = 1 the Lerch transcendent reduces to 

( 1 « , 

(8) ^ . 2Fl ( l + « ^ 

$(z, l,v 



and when n € {0, 1,2,... .}, $(z, n, v) has the hypergeometric representation |19) 

Vn 

1 +v n 



(9) *(z,n,v) = v-^F n ( 1 ^ - 



yielding 

(10) C(n,«) ^i^f 1 



and 



(ii) 



\ / n+1 



1 l/2„ 



— n+1 



lri+1 

? 

— 7; 



and thus due to (JTJ and Q we have the hypergeometric transformation 
(12) n+1 F n ( 1 ^ ^ = (1 - 2"") n+1 F n ^ 



3/2 r , 



2„ 



where the argument absent in p F q is assumed to be 1 and the symbol c„ denotes a parameter vector 
of length n where each element is equal to c (e.g. 53 = [5, 5, 5]). 

1.1.3. The Hypergeometric Polylogarithm. The polylogarithm, also known as Jonquiere's function, 
is defined V|t| <l,n€ {0, 1,2,.. .} by 

= n+l F n ( J l+1 \t ) t 



. 2„ 

(13) 



fei njuw* 



t Z^fc=0 t {2)1 

+ X^°° .fc r(fc+i) 
t Z^fc=o t r(fc+2) 

b l^k=a (l+fe)- 
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The hypergeometric representation (|1 of Li n (t) = n +i F n [ * ' " "+ 1 |£ ] t = LifYt) where 

\ Ol . . . O n J 

a\ . . . a„+i = ln+i and b\ . . .b n — 2 n is nearly-poised of the first kind 41, 2.1.1] since a\ + b\ — 
3 = . . . = a n + b n = 3. The notation Li n (i) refers specifically the hypergeometric form of Li n (t). 
The derivatives and integrals of Li n (i) satisfy the recurrence relations 



/o^dfl ^ =Li n+ i(*) 

(15) »+^« ( I s ) d * - n +2 F n+1 ( ^ |t) * 

and the reflection functional equation for Li„(l) = is 

(16) Lin(l) =^fefe% 

Li^(t) is seen to be (n — l)-balanced (|117p with the trivial calculation 

n n+1 

(17) ^2-^l = 2n-(n + l)=n-l 

k=i fe=i 

The usual defintion of Li n (£) requires analytic continuation at t — 1 but this is not necessary because 
the hypergeometric function converges absolutely on the unit circle when it is at least 1-balanced 
(fTTTjl which is true Vn ^ 2. The only Saalschiitzian polylogarithm is Li 2 (t) [32l Eq3.8] [20l 25:12] (2(3 
1.4.2] 

1.1.4. TTie Differential Equation Solved by Li^(i) and Some Combinatorics. Some combinatorial 
functions need to be defined before writing the differential equation solved by Li^ (t) . Let a partition 
be an arrangement of the set of elements 1, . . . , k into n subsets where each element is placed into 
exactly one set. The number of partitions of the set 1, . . . , k into n subsets is given by the Stirling 
numbers of the second kind 2, 1.1.3] 38, 2.7] defined by 

k \ _ „ n j k ( n 
n } - ^ j 

I K . - V" 3"(-l)"~ J ' 

\ 10 ) - l^j=o r(j+i)r(n-j+i) 

The kFk-i representation of / ^ 1 is (n— fc)-balanced (111 7[) since (fc— 1)— ((1— n)+2(fc— 1)) = ri— fc. 
The r-restricted Stirling numbers of the second kind < ^ I , or simply the r-Stirling numbers, 



n 

\. j r 

counts the number of partitions of the set l,...,n into k subsets with the restriction that the numbers 
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1, . . . , r belong to distinct subsets. [33] The recursion satisfied by j ^ I is given by 



(19) 



k < r 



n 



( k-1 








\ n 




In — 


1} 



n > r 



{[ 77 = 777 
77, 7^ m * S ^ 6 Kronecker delta. Specifically, the 2-restricted Stirling num- 

&ers [T51 A143494] appearing in the differential equation for Li^(rj) are given by 



r k i 




r fc i 


hi 


r jb-i i 


l « J 


L H 


l « J 




l « / 



\n-2 



(20) =^ES(-iy- fc (*7 2 )c? 

- f-1 ^ V fc ~ 2 (3+2)"- 2 (-l)^ 

- V 1 J 2-rj=0 j!(fe_2-j)! 

The (n + l)-th order hypergeometric differential equation (|119[) satisfied by f(t)=Li^(i) ([1 

r /(*) + ^/c*)(* 2 - *) ^=o 

(21) = \ m + (^/w) (V*- 1 { :+ 1 } - * m - 2 { n m + 1 }J ^ i 

has a most general solution of the form 

(22) f(t)=x + yG n (i) + Er=\ *m ln(i) m 

where x, y, Zi, . . . , z n —\ are arbitrary parameters and G n (t) satifies the recursion 

(23) G n (t) = 



ln(t - 1) 77 = 1 

Li 2 (l - t) + La(t - 1) ]xi(t) 77 = 2 

/ ^f^d< 77 ^ 3 



which has the explicit solution 



t 



i-t n ^ 

ln(t-l) n = l 

( 24 ) G «W = ) Li 2 (l -t) + ln(t- l)ln(t) n = 2 

( l n (t-l)-l n (l_t))l n (t)"-i ^ m(t)"- 2 t ■ m „ 

T(n) 6 T(n-l) " ? z 

The indicial equation of (j2Tj) at the t = 1 is 

(25) ind(LC(t)) = _ t(-i)-- 1 r(n-i-t)( t -n + i)» 

The (tt. + 1) roots of ind(Li^(£)) are the exponents of (|2ip which are simply 

(26) {t : ind(Li£(t)) = 0} = 0, 1, ...,77- 1, 77- 1 

where the last root n — 1 of ind(Li^(i)) is the balance of Li^(£) (IT7|) having multiplicity 2 thus 
inducing the logarithmic terms of (|22p . [16[ 15.31 and 16.33] These equations were derived by writing 
the differential equation for increasing values of n and then noticing that the developing pattern 
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of coefficients were combinatorial. After deriving the general combinatorial differential equation, it 
was solved for increasing values of n which resulted in nested integrals of prior solutions and then 
the general solution was derived from that pattern. 

1.1.5. The "Hypergeometric Form" of the Zeta Function. The main focus will be on the special 
case Li^(i) at unit argument where it coincides with the Riemann Zeta function at the integers. 
As with Li^(t), the symbol C F ( n ) refers specifically to the hypergeometric representation of £(n) 
at non-negative integer values of n. Using ([5]) and (fT3|). it can easily be seen that £(n) can be 
expressed as a generalized hypergeometric function (|116[) with 

( F (n) =L£(1) 
r ( 

-l-fn 



9 

- 2^ k =o fc] n?=i(2)* 

(27) 

~~ l^k=0 k\{2)l 

Eoo r(k+l) n 
k=0 r(fe+2)™ 

= C(M) 

= cw 

The value C F (0) =i F ( 1 |l) is singular and so must be calculated with the reflection equation 
([16]) to get Li^ (-1) =i F ( 1 | - l) = -| = C(0) which agrees with the integral form of ({t)Vt ^ 1 

(28) C(*)l t=0 = (| + th + 2 Jo°° sio(sarct ^ifff + * 2) ~ f ds) = Li (-1) = -| 



2. Number Theory, Continued Fractions, and Fractal Strings 

2.1. Fractal Strings and Dynamical Zeta Functions. A fractal string is defined as a nonempty 
open subset of the real line £1 C ]R which can be expressed as a disjoint union of open intervals Ij 
being the connected components offi. [231 3.1][50 1 [2l^ [ 2l2 1 [lT] [ 2"3 ] 

(29) n =U£i*J 

The length of the j-th interval Ij will be denoted lj . It will be assumed that f2 is standard, meaning 
that its length is finite, and that £j is a nonnegative monotically decreasing sequence. 

(m \n\ d = Z7=i^) d <^ 3d >v 

1 ' £ 1 >i 2 >...>i j > £ j+ i >--->o 

where 3d > means there is at least one value of d for which the statement is true. It can be the 
case that lj = for some j in which case lj is a finite sequence. The sequence of lengths of the 
components of the fractal string is denoted by 

(31) C ={ii)f =1 

The boundary of Q, in 1R will be denoted by dfl — K C which will also denote the boundary of C. 
Any totally disconnected bounded perfect subset K C 1R, or generally, any compact subset K C R, 
can be represented as a string of finite length |fi|x. A subset if of a topological space fi is said to 
be perfect if it is closed and each of its points is a limit point. Since here fl is a metric space and 
K C fi is closed, the Cantor-Bendixon lemma states that there exists a perfect set P C K such 
that K — P is a most countable. [35l 2.2 Exl7] As such, f2 can be defined as the complenent of 
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K in its closed convex hull, that is, — Q(K) is the smallest compact interval [a, b] containing K . 
The connected components of the bounded open set f2 = (a, b)\K are the intervals Ij of the fractal 
string C associated with K. 

2.2. The Gauss Map h(x). Let ti h = (0, l)\dfth where dti h = {±00, 0, ± : n e 7L\ is the set of 
discontinous boundary points of the Gauss map h(x) G O/iVx dClh, also known as the Gauss 
function or Gauss transformation, which maps unit intervals onto unit intervals and by iteration 
gives the continued fraction expansion of a real number 

I _ [II 

(32) ~7^T" 

= {x L } 
— - mod 1 

where [x\ is the floor function, the greatest integer^ x and {x} = x — [x\ is the fractional part of 
x. [H 2.1,3.9.1,9.1,9.3,9.7.1] '42, A.1.7] Clearly h{x) is also defined outside of fi 

1 



h(x) 



(33) 



h(x) 



h{x) 
i + 1 



x > 1 
-1 < a; < 1 
x < -1 



since 



(34) 




As can be seen in Figure [I] h(x) is discontinuous at a countably infinite set of points of Lebesgue 
measure zero on its boundary dQh 



(35) 



y : lim h(x) 7^ lim h(x) > =dQ 



x-yy 



x— }-y~< 



±00, 0, — : n e % 



The left and right limits of h{x) when x approaches an element on the boundary dflh is given by 
(36) 



2.2.1. The Frobenius- Perron Transfer Operator. The Frobenius-Perron transfer operator [421 2.4.4] [3T1 
2.3.3] |T3l 1.3.1,8.2] [40l 1.8,2.4] of a unit interval mapping f(y) describes how a probablility density 
p(y) transforms under the action of the map. 

(37) [U fP ](x) =jS(x-f(y))p(y)dy 
where 8 is the Dirac delta function and 6 is the Heaviside step function. 

J5(x)dx =6{x) 

(38) a . , _ / x < 



9(x) 



1 x > 



The function f(y) is the map being iterated and p{y) is some density on which the transfer operator 
U acts. Essentially, iteration of the map transforms points to points and iteration of the transfer 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



_x 

{h(x):x=Q.l..l} 



Figure 1. The Gauss Map 



operator maps point densities to point densities. The Gauss-Kuzmin-Wirsing(GKW) operator is 
obtained by applying the transfer operator to the Guass map. (44] 2] [50] [46] 

( 39 ) WW =E^r^# 

By changing the variables and order of integration in (|65[) an operator equation for £(s) is obtained. 
C(s) =^ l -s^x[U h x s - l ]dx 

= t=T ~ s Jo x 1 5 ( x ~ HvW^dydx 

(40) = - a £ x J S(x - (y-i - [y^\ ) )y^dydx 

S fl V^°° ( n + x ) J 

— ~ ~ S JO X 2^n=l (n+x) 2 aX 

An operator similiar to is 

(41) P*p](*) =En =1 P^)-p(^, 
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The action of [Uhp](x) on the identity function x — > x is given by 

[U h x)(x) = En=i\^+fy 

Eoo 1 

(42) = _£%ti2 

( 1 x+1 x+l x+1 
" : ' x + 2 x + 2 x + 2 



~~ (x+i) 3 

where ty( n '(x) is the polygamma function (|122l) . The area under the curve of [Uhx](x) is 

(43) [U h x](x)dx = L^i dx= _ 

The identity action of [Shp]{x) is 

(aa\ [Shx](x) = En=l n ~ iii 

1 ' = J + ^(x + l) 

where 7 is Euler's constant 

7 = lim^oo J2k=i \ - m ( n ) 
= \im s ^i C(s) - 

= lim^i C(s) +j f 1 ° h(x)x s ~' 1 dx 

= lim s ^i j^- — s h(x)x s ~ 1 dx + h(x)x s ~ 1 dx 

= 1 — J Q h(x)dx 

w 0.57721566490153286060651209 
and the area under its curve is given by 

(46) fi[S h x](x)dx = f£>y+9(x + l)dx=l-'y 

2.2.2. Piecewsise Integration of h(x). The Guass map h(x) £ O/j is piecewise monotone [101 2.1] 
between the points of dflh, and thus partitions the unit interval infinite covering set of decreasing 
open intervals seperated by d£lh- 13, 5.7.1] Let I n be an infinite set of open intervals 

(1, 00) ri = 



(47) In = 



< n < 00 



It is easy to see that 
(48) 



(O,O) = n = oo 



n h \jasi h = [0,1] = IX=i 7 « 

[0.OO] =IX=o'n 



Define the Gauss map partition h n (x) where {/i„(a;) 7^ : x G /„} as a piecewise step function 



(49) ' 



l-an 1 <X< 1 



otherwise 

/ 1— xn \ I n I xn+x — 1 \ n ( xn— \\\ 
[ { n+1 ) - J 
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where 6{t) is the Heaviside step function We can reassemble all of the {h n (x)}^ =1 to recover 
h(x) 

h(x) = J2n=i h n {x) 

(50) =E ^ i (l ¥B )( tf ( 2 ^i)_ tf ( S 

where only one of the h n (x) is nonzero for each x. By setting n = I ij in P^|) we get 

f ___LiJ i <r< i 

(51) [ otherwise 



UJ+w ~ \UY 

Define the partitioned integral operator [Pf(x);x] (n) by 

(52) [Pf(x);x](n) =[Pf(x)](n) = [P/] (n) = /(a;)ds = /(a:)da 

where by convention we have 

, , [Pf{x)](0) =]hn n ^ 0+ J I J(x)dx =f™f(x)dx 

[Pf(x)](oc) =lim n ^ 00 J I J(x)dx =J°f(x)dx =0 

Thus 

( c 4) /o/tod* =E~i[i'/W](n) = E£= 1 /, B /(*)<!* 

1 j r/W^ =Zn= lPf(z)](n) =E^ol/(^ 



Each interval /„ has the length 



£J„ =[Pl](n) 
= f\_ Ida; 



(55) = 

n ^ n+1 
— n(n+l) 

The elements n(n + 1) are known as the oblong numbers [T5, A002378]. It is seen, together with 
(PI) , that 

|fi h | = Jo Ida; 

— Z^n=l "n 

Eoo 1 
71=1 



(56) ^n=l n(n+l) 

K ' =\Shx]{l) 

= 7 + ^(2) 
= 1 

The piecewise integral operator [Pf(x);x] (n) can be used to calculate the area under the curve 
of h(x) which is also equal to the area under the curve of [Shx](x). Let the length of the ra-th 
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component h n (x) be denoted by 

lh n = [Ph(x);x](n) 
= Jj h(x)dx 

(57) = //" K(x)dx 

= Jo h n(x)dx 

_ ln(n+l)n+ln(n+l)-ln(n)re-ln(n)-l 
— n+1 

Regarding h(x) as a fractal string Ch = {h n (x)}^ =1 its length \Ch\ is given by 

\£h\ =f h[x)dx 

— En=l 

ln(n+l)n+ln(n+l)— ln(n)n— ln(n)— 1 

(58) " tt n=1 ™+i 

= Jo [S h x]{x)dx 

= f Q 1 'Y + y(x + l)dx 

= 1-7 

If n = in (|48p we get the interval 1$ = (y, ¥) = (1, oo) and 

(59) £/io = J Io h(x)dx = /~ idx = 
but if we choose a finite cutoff then 



ln(y) 



(60) 
and 

(61) r = j 
v ; Jl y « » ( n -l) 2 

thus 

n=2 Jl d V =J2n=2 Jl 1± ^ d V 

n=2 Jl -ir d V 

(62) - V°° 1 

= C(2) 

= 

6 

2.2.3. T/ie Meffin Transform. The Mellin transform [Ml H.10.8]|3| 3.6] is defined as 

(63) Mh b }f(x) =J b a f(x)x s - 1 dx 

where the usual definition of the Mellin transform is Mx^f^f(x). Somewhat incredibly, by taking 
the Mellin transformation of h(x) over the unit interval, we get an analytic continuation of £(s) 
which is convergent when s is not equal to a negative integer, 0, or 1. When s is a negative integer 
or the limit or analytic continuation must be taken since the series is formally divergent at these 
points, and of course the series s = 1 diverges. [25] [H] [4"5] 

M^ s h{x) = [Ph(x)x s - 1 ;x] (n) 

(64) =fi T (i-H\)x s - 1 dx 

_ n(n+l)~ B +s(n+l)~ B -n 1 ~ B 
- s(s-l) 
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= — s J Q h(x)x s ~ 1 dx 
(65) =^-'/o 1 (§-|iJ)* , - 1 d* 

s Y^oo n(n+l) s +s(n+l) s — n 1 3 

~~ I^T _ S 2^n=l s(s-l) 

The term -^j changes to ^ = - (-1) by subtracting the residue 10.41] [48j 6.1] of 

1 



(66) M^% s h(x) = / h(x)x s - 1 dx = / dx 



s- 1 

1 

at the singular point s = 1, which happens to coincide with — 



(67) 



Z^s=2 
-1 



2.3. The Harmonic Sawtooth w(x). Define the harmonic sawtooth map w(x) £ ilh\dflh which 
shares the same domain and boundary as the Gauss map h(x) to which it is similiar, and also 
has the property that its Mellin transform is the (appropriately scaled) zeta function. The n-th 
component w n {x) is defined over the n-th interval 



(68) Wnix) 



n(xn + x — 1) — 7T < x < — 

\ / n+l n 

otherwise 



= n (xn + x-i)(e (^fr 1 ) - & (^)) 

and by the substitution n — > [-J we have 



w { x ) =Y,n=l W n(x) 

(69) = J2™=i n ( xn + x - 1) 

= L^J (*LSJ+*-i) 



n+l 



Unlike which is nonzero outside of |x| > 1, the (harmonic) sawtooth map has w(x) = 0V|x| > 1. 
The length of each component of w(x) is 

£w n = [Pw(x);x] (n) 
(70) = J In w(x)dx 



2(n+l)n 

So that the total length of the harmonic sawtooth string C w is 

|Au| = J 1 w(x)da; 

- V 00 £w 

(71) _4<Sg =1 i 

2^n=l 2(n+l)n 
~~ 2 
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X 

{w(x):x=Q.l..l} 



Figure 2. The Harmonic Sawtooth 



The infinite set of Mellin transforms of w n (x) 



(72) 



= [Pw(x)x s 1 ;x\ (n) 
= f " n(xn + x — l)x s dx 

n+1 

= J* n(xn + x - V ' a I £S±x-l 



n+1 



i(n+l)~ a +s(re+l)~ s -ra 1 
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are summed to get an expression for {((s) : 9t(s) ^ N - } 

(73) = En=i s 7 ± jI^ fT n ( xn + x - ^x^dx 

Eoo / n(n-\-l)~ s -\-s(n-\-l)~ s —n 1 ~ s 

n=l 6 s-1 ^ s(s-l) 

Eoo n(n+l) -s — n 1 ~ s +sn _s 
n—1 s—1 

2.4. The Prime Numbers. Let P = {2,3,5,7,11,13,17,19,23,29,...} denote the set of prime 
numbers and Ni = {1,2,...} and N = {0, 1,2,.. .}, N = {...,-2,-1,0,1,2,...} be the set of 
positive, non-negative, and signed integers. 

2.4.1. The Prime Counting Function: tt(x). The prime counting function tt(x) counts the number 
of primes less than a given number. It can written as 

(74) *(x) =ESJl 

which is essentially a step function which increases by 1 for each prime. [9l 15.11] 

2.4.2. von Mangoldt and Chebyshev's Functions: A(x),0(x),tp(x). Chebyshev's function of the first 
kind 0(x) is the sum of the logarithm of all primes ^ x 

6{x) =Y.tlHp k ) 

(75 

where pk € P is the fc-th prime. [71 4.4] The generalization of n(x) is the Chebyshev function of 
the second kind 

=E^' r£Nl} ln(p) 

= Efa 2(x)J *(**) 

(76) = ln(lcm(l,2,3,..., [x\)) 

= x _ Ml^!) _ ln(27r) _ E C(p)=o ^ VI(p) ^ Q 

where the first sum ranges over the primes p £ P and positive integers r and the sum over p is 
von Mangoldt's formula where p ranges over the non-trivial roots of £(s) in increasing order. The 
function lcm(. . . .) is the least common multiple, and A(x) is the von Mangoldt function. 

ln(p) {n — p k : p £ P, k G Ni} 

otherwise 

/ lcm(l,2,...,n) 

1 lcm(l,2,...,n-l) 



A(x) = 

(77) 



In 



A(s) is related to ((s) by 

(78) - = V ^^Vft(s) > 1 



^C(s)_y>AK 



n=l 

Chebyshev proved that tt(x), 9{x), and VK 2 -) have the same scaled asymptotic limit. 
(79) hm . v \ = hm — — = hm = 1 



X— >OQ 



(ln(z)) 
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QH 1.3] QZ1 1.4] 4.3-4.4&3.1-3.2] [5] 15.11]. Note that [5] incorrectly defines V(ar) as 
ln(gcd(...)). 

3. Analytic Continuation 

3.1. Continuation of n +iF n Near Unit Argument. The continuation formula for Gauss's hy- 
pergeometric function 2F1 near unit argument is well known 

r(ai)r(q 2 ) p / ai a 2 , \ _ ^00 (_l)"(l- z )" r(ai+n)r(a 2 +n)rQi-rt) 

(80) r(bl) 2 1 ^ 61 1 y ~ ^"=0 ™ ! r(a 1+Sl )r(a 2+Sl ) 

1 /-i _ \si r(ai+si+w)r(a 2 +si+Ti)r(-si-n) 
>\ L z ) Z^n=0 n! r(ai+si)r(a 2 +si) 

where s± — b± — a\ — ai is the balance (|117l) of 2-F1 which must not be equal to an integer, that 
is, iF\ cannot be si-balanced. A function is said to be fc-balanced only when k is an integer. When 
IZ(si) > the value at z = 1 is finite and given by the Gaussian summation formula 

, Fl ( ai ° 2 

(81) V bl / _ rC^-oi-oa) 

r(oi+si)r(o 2 +si) 

It is obvious that limt^i Lif (t) — Umt^i 2-P1 ^2 1^ = £ F (-0 — 00 is 0-balanced and of course 

equal to the divergent harmonic series so the continuation formula does not apply. However, Biihring 
and Srivastava [6] [5] generalized this relation to all n +iF n by expanding (|8Tj) as a series then 
interchanging the order of summations to derive a recurrence with respect to n 



(82) 



p / oi,...,a n+1 I \ _ r(6 n )r(&„_!) 
n+1 71 V 6i,...,i>„ ' y ~~ r(a„ +1 )r(b„+6„_ 1 -a„ +1 ) 

(b,i-a„+i) m (b,i-i-an-i) m p / dl, . . . , a n , 
' Lm=0 ( 6 „ + b „_ 1 -a„ +1 ) m m! { b ly . . . , 6 n _ 2 , &«-! + &„ - «„+! + m ' , 

which is valid \/{lZ(a,i) > : 1 ^ i ^ n + 1}. The m-th term of the summand in (|82p is contiguous 
f|4. 1 . 2[) to the (m — l)-th and (m + l)-th terms and thus a linear relationship can always be found 
between neighboring terms. 

3.2. The Continuation of Li^ (t) and £ F (n) via Contiguous Functions. There are 4 functions 
contiguous (|4.1.2[) to hi n (t), only 3 of them are unique, and just 1 of them is interesting. The 
functions are obtained by shifting one of the numerator parameters a* ± 1 or shifting one of the 
denominator parameters hi ± 1. Shifting any of the a parameters or any of the b parameters will 
suffice since they are all equal and p F q is invariant with respect to the ordering of parameters. Let 
and denote the parameter vector c n where one element is shifted up or down by 1. 

c+ = c„_i, c +1 = c, . . . , c, c + 1 

(83) _ 1 1 

c n = c n _i,c- 1 = c, ...,c,c- 1 



n-l 
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For example, 4g~ = 4,4,5. Two of the four functions contiguous to Li^(t) are identical 
(84) n+iF n ( \t) t = n+1 F n ( |A t = Ld(t) 



Shifting any <Zj up is equivalent to shifting any hi down, both operations take Li^(i) to Li n _ 1 (f). 
Shifting any down results in the identity function since it puts a in the numerator. 



(85) 



n+l-^n 



lfi+1 



+ 1 It t = t 



Thus, the only interesting function continguous to Li^(i) is obtained by shifting one of the denom- 
inator parameters up. Let this function be denoted by Li^ +1 (t) 



(86) Li£ +i (t) = n+1 F n ( t^ 1 I-/ ! = 



h (2Vt) - j-h {2V~t) n = 

T-T- 1 "=1 
(_!)« (i _ + Li fc (t) ) ? i > 2 



where I n (x) is a modified Bessel function of the first kind |34j 65] jlOi 6.9.1] 



(87) 



In{x) 



r(n+l)2" F l I n + ! I 4 



Before applying (j8"2")l . the notation will be simplified by extending (|83l) so that repeated shifts can 
be written more easily 

& =Zn-l,C + j = C,...,C,C + j 

(88) -t_j _ . _ 

c n — c n _i, c j — c, . . . ,c , c j 

71-1 

where clearly = c^" 1 and = c^ 1 . The goal is to extend Li^ +1 (t) to all Li^ +m (t) by repeated 
application of c^ 1 . Applying (j52"j) to ([TB")) gives the continuation of Li n (t) — > Li^ +1 (t)Vn ^ 1 by 
setting ai... n +i = l n +i and = 2„ which results in 



(89) 



- tV 00 



tV°° 
Z^m=0 



i-n+1 



|t t Vn ^ 



1, 



gn-2, 3 + m 

(m+l)(m+2) 

V 

(m+l)(m+2) 



Vn > 2 



Vn > 2 



since 



(90) 



r(6 n )r(6„_ 1 ) 



r(2)r(2) 



r(a„+i)r(b„ + b n -i - a„+i) r(i)r(2 + 2-i) 2 



TWO NEW ZETA CONSTANTS: FRACTAL STRING, CONTINUED FRACTION, AND HYPERGEOMETRIC ASPECTS OF THE RIEMANN ZE1 



and 



(91) 



(b n — a n +l)m(b n -l — 0>n-l)7 



(l)m(l)r 



{b n + b n -x-a n+ i)m\ (2 + 2 - l) m m\ (m+l)(m + 2) 

The denominator parameters 2 J ^\ m in (|89|) are simply 



(92) 



n + l + m 
z n-l 



The numbers (to + 1)(to + 2) are known as the oblong numbers, [151 A002378]. By simply setting 
t = 1 in (|59")l we get the continuation from C^v?" 1 ) — > ( F ( n + l)Vn > 1 



/ 



(93) 



C F (ra) = n +l F n 



1)H 



m=0 



In )\ 



(m+l)(m+2) 



V 



Vn ^ 2 



/ 



The justification in saying that Li F (i) and C F ( n ) are continued to Li F +1 (i) and ( F (n + 1) comes 
from the fact that the first term in the summand of the continuation (f89|) from Li F _ x (t) — > Li F (t) is 

contiguous to Li F _ x (t), that is, „F n _i ^ i"j |ij is contiguous to Li^l^t) = „F„_i ^ i n |i 

The continuation formula f|89|) gives interesting answers for n = and n = 1 which suggest an 
alternative to "the analytic continuation" of ( (t) which is different from the usual 1 _ 1 2 -t S^Lo(^ n + 
We have 



/ 



C F (0) =E7n=0 



(m+l)(m+2) 



V 



(94) 



y^oo -(J m + i(2)m+/m+iC2)-/ m (2))r(m+3) 
2^m=0 (m+l)(m+2) 

/o(2) - 1 
1.2795853023360 



C F (i) =E™ =0 



\ 



(m+l)(m+2) 



V 



v-^co e(r(m+3)-r(m+2,l)m-2r(m+2,l)) 
2-~im=0 (m+l)(m+2) 

Ei(l)- 7 
1.3179021514544 



where Ei(x) is the exponential integral [lOi 6.9.2] 

Ei(.x) - ~ - ->- ^ 



(95) 



' 2^2^ Z^fc=i fer(fe+i) 



_ Info- 1 ) Info) 
' 2 ' 2 



x 2 F 2 



1 1 

2 2 
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and r (a, z) is the incomplete Gamma function 



(96) 



r (a, z) — r (z) 



a + 1 



So we have the "hypergeometrically continued" values C F (0) = ^o(2) — 1 and C, F {\) = Ei(l) — 7 
whereas the "real" values are £ (0) = — ^ and £ (1) = 00. In terms of reciprocal probability we have 



(97) 



C F (Or 1 S 78.15% 
C F (1) _1 = 75.88% 



3.2.1. Lif (t) Li F (i) and C[(t) Cf (*)• Tne continuation C F (t) from re = 1 -+ 2 via is 
straightforward 



(98) 



C F (2) =3^ 2 



1 1 1 

2 2 



Eoo 
m=0 



Eoo 
m=0 



2 Ft 



1 L ^ 

6 + m 



(m+l)(m+2) 



Y^cxj r(m + 3)r(fc+i) 

Z^k = Q r(m + fc + 3) 

(m+l)(m+2) 
1 (m+2) 



Eoo 
m=0 (m+l)(m+2) (m + 1) 

Eoo 1 
m=0 



(m+1) 2 



The continuation of Lif (i) to Li F (i) is a bit more complicated 

. ,-. /Ill 



Li£ (t) =3 ^2 



(99) 



then r 2(171, t) is given by 



2 2 



EOO 
m=0 



m + o 



(m+1) (m+2) 



V 



(100) 



r 2 {m,t) 



,p m ( _l ) " + l r(m + 2 )(^(, T1 -„ + l)-^(„ 1 + 2))(-l)'"e^("' + 2 )t" 
^" = r(n + 2)r(m-n + l) 



(t-iy 



(m+l)e*< m + 2 )t m + 1 
■ ln(l-t) 



m+1 



so Li F (i) is equal to 
(101) 



oo v^" 1 r(m+2)(^(m-n+l)-^(m+2))(-l)" 
2^n=0 r(n+2)r(m-n+l) 



m=0 



(m + l)e^( m +2)t m + 1 



(t- l) m+1 t- 2 - m ln(l-*) 

TO + 1 
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where %jj(m) = ln(lcm(l, 2,3,..., m)) is Chebyshev's function of the 2nd kind (|75|) and W(m) is the 
digamma function 



(102) 



^r(x) = — ln(r(x)) 
dx 



r(x) 



3.2.2. C F (2) C F (3). The continuation from ( F (2) to C F (3) via ([S3) is carried out like so 



C F (3) 



(103) 



m=0 



1111 

2 2 2 



/ I 1 1 1 

2 3 + m 



(m+l)(m+2) 



V 



C(2) 



/ 



Each term in the summand r 3 (m) has the form + 53(771) where of course £(2) = ^ and g 3 (rr 
is a rational function of to which follows a 3rd order linear recurrence equation [28j 8.2] given by 

(104) 



q 3 (m) = q 3 (to + 1) (to 3 + 8 to 2 + 21 to + 18) 

+<? 3 (m + 2) f-2 m 3 - 20 m 2 - 67 771 - 75) 
+<?3 (777 + 3) (to 3 + 12 to 2 + 48 m + 64) 



(105) <7 3 (m) 
The solution to which is given by 



1 m = 

I m = l 

■42- m — 2 

108 m ~ A 



(106) 

so the summand r 3 {rn) is 
(107) 

Thus (jllip is also equal to 

(108) 

Thus 



93(777) = m+1 



7 *3( m ) =lrfi+93( 



rn 



^ (1) (2+m) 
m+1 



c F o) =e:.„ 



00 ^W(m+2) 



m+1 



r3(m) 



_ ^< 1 '(m+2) 

m+1 
_ C(2,m+2) 
m+1 

— 6^ l^k=\ k 2 



(109) 



>J 7 2 



m+1 

1 to + 2 m+2 

777 + 3 771 + 3 



(m+l)(m+2) 2 
1 1 1 

2 3 + 777 

(m + l)(m+2) 
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The first 10 terms of {q^ (m) : m = . 

5 49 205 5269 

' 8 ' ~ 108 ' ~ 576 ' ~ 18000 ' ~ 21600 ' ~ 1234800 ' ~ 5644800 ' ~ 57153600 ' ~ 12700800 



(110) 



1. 



. 9} are 
5369 



266681 1077749 9778141 



1968329 



The denominator of (|110j) appears to be [HI Al 19936], the least common multiple of denominators 
of the rows of a certain triangle of rationals and the numerators are |15[ A007406], the numerator 
of E&=i W lrom (|123p which, according to a theorem Wolstenholme, p divides numer(<? 3 (p — 1)) 
where p € P is prime. [12] [4] [T] 

3.2.3. ( F (3) -> C F (4)- The continuation from ( F (3) to C F (4) via (|33) is given by 



C F (4) =5^4 



(111) 



1111 

2 2 2 2 
(.(I I 



E 



m=0 



1 

2 2 



(m+l)(m+2) 



V 



Em=0 r 4(m) 



The summand r4(m) has the form 
(112) r 4 (m) =a(t,m) 



b(t, m) 



H(m+1) Li 2 (t) 
(m+l)t 



Li 3 (t) 



r (m+l)t 

where a(t, m) is an (m + l)-th degree polynomial and b(t, m) is a (m + 2)-th degree polynomial(thc 
determination of which is left to an excercise for the reader or the topic of another article, but is 
readily obtained with the help of Maple [2"7]). and H(n) is the n-th Harmonic number 

H(n) =E? =1 ± 

= + 1) +7 

(113) = Efc=l fe2+fc„ 

_22_ P | 1 1 
«+l3 2 \ 2 



n 
n 



The polynomial 6(t, m) vanishes when i 

C(4) 



1. An interesting set of formulas for £(4) is 



_ v^°° 1 ,(2) (rt+l)+2C(3) 
~ 2^n=l 2n(n+l) 



(114) 



n=l 



2n(n+l) 



4. Appendix 

4.1. The Generalized Hypergeometric Function : p fg. The Pochhammer symbol is defined according 
to 



(115) (n) fc 

The generalized hypergeometric function [39] [ 
finite products of Pochhammer symbols 



T(n+fc) 
: T(n) 

4.1] is defined as an infinite sum of quotients of 
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The function p F q is said to be fc-balanced [5] if the sum of the denominator parameters b\ . . .b p 
minus the sum of the numerator parameters a\ . . . a p+ i is an integer. 



(117) k = ba\( p F q ) = 5 « - 



q 

n— 1 n— 1 



The value k is the characteristic exponent of the hypergeometric differential equation at unit argu- 
ment which is equal to the maximum root of the corresponding indicial equation and so determines 
the behaviour of the function near this point. A 1-balanced function is said to be Saalschiitzian. 
[33 2.1.1] 



4.1.1. The Differential Equation and Convergence. The function p F q converges when 
(118) 



p ^ q V|t| 7^ oo 

p = q + i y\t\<i 

{p = q+l :bal( p F 9 ) > 1} V|i| = 1 

p > q + 1 yt = 

where bal( p F q ) = X)«=i &n — ELi fl « * s * ne balance of the parameters (|117p . The differential 
equation solved by p F q is of order max(p,q+l) 

(119) UtHiOt + b.-^-tHiet + a,) ] /(t) = 

where /(i) = p F g ( , p |t ) and 6* t = t^ is the differential operator. When p = q + 1 

\ "i , . . . , o q J 

(Till)]) has the form 

(120) oo/(t) + t 9 ^+lf(t) + E * n_1 ( ta » - 5 ")dF / W = 

n— 1 

Hi 4.2] [211 Ch3][3j 44-46] [4TJ 2.1.2] 

4.1.2. Contiguous Functions and Linear Relations. Any two hypergeometric functions p F q {a..., 6...; z) 
and p F q (c..., d...; z) are said to be contiguous if allp+g pairs of parameters (ai, ci), . . . , (a p , Cp), di), . . . , (b g , d q 
are equal except for one pair which differs only by 1 . There are 2p + q linearly independent rela- 
tions between the (2p + 2q) functions contiguous to p F q (a,_,b,_; z) where the relations are linear 
functions of z and polynomial functions of the parameters a... When any {a^ = dj : i =/= j} 
or {bi = bj : i ^ j} in p F q (a..., z) there will fewer unique contiguous functions than if all the 
parameters were unique since the hypergeometric function is invariant with respect to the ordering 
of parameters. gD 2.2.1] [31 48] gH] QUI 4.3] gH] [53] 

4.2. Other Special Functions. 

4.2.1. Polygamma ^^ n \x). The polygamma function is the n-th derivative of the digamma (|102p 
function 

(121) = ^*(z) 
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*<")(.t) = 

(122) 



and is defined as an infinite sum, a Hurwitz zeta function (U), or a hypergeometric function when 
x is positive integer 

(E£U I - fe+VT^r) - 7 » = 
V°° " ! (-i) n n > 1 

^3«( 1 1 :+i))-> 

( n = 

~ I (-l) Il+1 n!C(x,n + l) 1 

or as a finite sum when x is a positive integer and n = 1 10, 1.16] 

(123) = if - Efci w 
4.3. Notation. 

Z {...,-2,-1,0,1,2,...} 

N {1,2,3,....} 

(124) Nx- {...,-3,-2,-1} 

N {0,1,2,3,....} 

N - {...,-3,-2,-1,0} 
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